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Abstract. Genetic algorithms, as implemented in optimal control strategies, 
are currently successfully exploited in a wide range of problems in molecular 
physics. In this context, laser control of molecular alignment and orienta- 
tion remains a very promising issue with challenging applications extending 
from chemical reactivity to nanoscale design. We emphasize the complemen- 
tarity between basic quantum mechanisms monitoring alignment /orientation 
processes and optimal control scenarios. More explicitly, if on one hand we can 
help the optimal control scheme to take advantage of such mechanisms by ap- 
propriately building the targets and delineating the parameter sampling space, 
on the other hand we expect to learn, from optimal control results, some robust 
and physically sound dynamical mechanisms. One of the basic mechanisms for 
alignment (i.e., molecular axis parallel to field polarization) is related to the 
pendular states accommodated by the molecule-plus-field effective potential. 
The laser control of alignment can be reached by an adiabatic transport of an 
initial isotropic rotational state on some pendular state trapping the molecule 
in well-aligned geometries. Symmetry breaking mechanisms are to be looked 
for when orientation (i.e., molecular axis in the same direction as field polar- 
ization) is the goal of laser control. Two mechanisms are considered. The first 
is based upon an asymmetric pulse combining a frequency u> and its second 
harmonic 2uj resonant with a vibrational transition. A much more efficient 
mechanism is the so-called "kick" that a highly asymmetric sudden pulse can 
impart to the molecule. Half-cycle pulses, within the reach of current experi- 
mental technology, are among good candidates for producing such kicks. Very 
interestingly, an optimal control scheme for orientation, based on genetic algo- 
rithms, also leads to a sudden pulsed field bearing the characteristic features of 
the kick mechanism. Optimal pulse shaping for very efficient and long-lasting 
orientation, together with robustness with respect to temperature effects, are 
among our future prospects. 



1. Introduction 

Laser-induced molecular alignment and orientation are challenging control is- 
sues with a wide range of applications, extending from chemical reactivity to 
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nanoscale design ^ |3 I20L I24L I26L 12 7| . They address molecular manipulation, 
involving external angular degrees of freedom, aiming at a parallel positioning of 
the molecular axis with respect to the laser polarization vector (alignment) or, 
even more demanding, with a given direction (orientation). One of the basic mo- 
tivations remains the drastic increase of reactive cross sections, highly sensitive to 
stereodynamical effects (frontal collisions). Laser-induced isomerization jlOL I19j . 
isotope separation molecular trapping |23j . or surface processing and catal- 
ysis [3 I24L 12 7j are some of the illustrations of the widely growing interest in 
molecular manipulation. 

The most common idea is to excite a broad rotational band A J in order to re- 
cover a narrow angular distribution A^, through the celebrated relation A JA6' ^ h 
(J being the total angular momentum and 9 the molecule-laser angle). Our ap- 
proach to this problem is, however, different. The emphasis is put on the depiction 
and study of some laser-induced dynamical mechanisms that are relevant to the 
alignment/orientation process: pendular states in high-frequency fields; two-photon 
(fundamental and second- harmonic) excitations; the kick mechanism. The question 
of how to take advantage of such mechanisms in a control scheme is brought up only 
in a second step. Adiabatic transport from an isotropic distribution to a pendular 
state; two IR laser pulses with frequencies in a ratio 2, the second harmonic be- 
ing in resonance with a vibrational transition; very asymmetric and short-duration 
(sudden) half-cycle pulses imparting sudden kicks to the molecule, are some of the 
scenarios which are depicted and discussed in the following. Finally, in a third 
step, this knowledge is fully exploited in optimal control schemes, not merely for 
initial guesses in parameter sampling, but mainly to get a better understanding 
and interpretation of otherwise black-box type calculations. 

The paper is organized as follows. Section |2 is devoted to intense laser align- 
ment of the HON molecule, through pendular states as the elementary mechanism 
and adiabatic/sudden transport as the control strategy. Section |3 presents the two- 
color orientation of HON by a parity mixing scenario in the rotational distribution. 
The kick mechanism with a half-cycle pulse sudden excitation scenario is examined 
in section 01 as illustrated on the LiCl molecule. Sudden-excitation-based opti- 
mal control scenarios, with a thorough analysis of the role of different orientation 
criteria, are analyzed, using HON and LiF, in sectional 



2. Laser-controlled alignment 

2.1. The model. The linear HON molecule, taken as an illustrative example, 
is described by its internal stretching motions (collectively labelled R) and external 
rotational motions described by its polar 9 and azimuthal tp angles with respect to 
the laser polarization vector (see figure^. The Hamiltonian involves a single Born- 
Oppenheimer potential curve VCR) (the molecular electronic ground state) and the 
radiative coupling _ffi.ad(R', ^, t) given by the scalar product of the field- induced 
dipole moment /.t(R) with the electric field amplitude £(t), i.e., 

(2.1) HCR,9,^,t) =fR+ ^^J^ + VCR) + H,.,^(R,9,v,t), 

where Tr is the kinetic energy operator in R, the angular momentum operator, 
and /(R) the moment of inertia of the molecule. The dipole moment itself is taken 
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Figure 1. The linear HCN molecule positioned, with respect to 
the polarization axis £ of the laser, by its polar 9 and azimuthal 
angles. Internal stretching coordinates are collectively noted R — 
iR,r). 



as the sum of a permanent dipole /Xq and a polarizability a 



(2.2) 



such that 



(2.3) 



IJ.{R,t) = fi„ + ^a -Sit), 



1 



i?rad(R,6',<y5,i) = -^^a£{t) cose--£\t) (a|| cos^ 



a I sm 



In the last equation, we have made use of the fact that, for a linear molecule, 
a = 'D\a.g{axx,ctyy,azz\, with axx = otyy = ai_ and azz = Oi\\, when the z axis is 
chosen to be coUinear with the molecular axis. Due to cylindrical symmetry and 
since the laser-molecule coupling term does not depend explicitly on tp, the motion 
on if can be considered separately and is no longer taken into account. 

The laser pulses are of the form £{t) = £o{t) cos {uit + (j)), with £o{t) the time 
varying envelope of the electric field of frequency u) and phase (j). We will consider 
laser pulses of a maximum intensity of the order of 10^'^ W/cm^, to avoid ionization 
damage to the molecule. The frequency uj is taken in the infra-red region, so it is 
considerably higher than the rotational frequency of the molecule. 

The dynamics of the molecular system are then obtained by solving numerically 
the time-dependent Schrodinger equation 



d 

ih—ip{R, 0, t) = H(R, 0, t)iPCR, 0, t) 



(2.4) 

(on a discretized grid using the split-operator method |3l ll3L[l4] ). starting from an 
initial wave function ip{R, 6,t = Q) taken, in most cases, as the ground (rotationally 
isotropic) state of the molecule. 

2.2. Pendular states. In this high-frequency regime, two different motions 
are considered: a fast motion associated with the electric field wt and a slow motion 
associated with the molecular rotation described by 0. To simplify the presentation, 
we freeze the molecular internal motions, assuming a rigid rotor behavior. The 
eigenvectors of the slow motion £,nj{0) obey, after averaging over the fast motion 
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within an optical period (and some high-frequency approximations) a Schrodinger 
equation given by jl8| 

!i2 _ t,2 n,f2 / ,,2 \ p„ A^,2f-4 



— H ^ + Aa + — sin^ 6* H 2^ sin^ 26* 



2/ " '2/ sin^ 61 V /c^ V 4 256/tj2 

(2.5) ^ (a„- Jfiic> + ^ )e„j(^), 



where Aa = a|| — a^. The ^-dependent laser-induced effective potential has a sin^ 9 
shape and supports the eigenvectors £,nj{0) which are localized around 6 = Q and tt 
(aligned geometry), in a way which is more marked for increasing intensity. Because 
of the resulting libration motion for a molecule that is in one of these eigenstates, 
they are dubbed pendular states |15L 12 9| . 

It is interesting to note a close analogy with the classical Lagrange equation of 
the laser-driven rigid rotor under the combined effect of the permanent dipole and 
the polarizability |12| . 

. ^ cose / f^n \ £n . Aa^f^ 

2.6) 76' ~— = Aa-f^ ^sin26>- sin46l. 

The classical force of the right-hand-side of equation H2.6f) is nothing but the gra- 
dient of the effective potential of the Hamiltonian version, equation H2.5|l . This 
classical analogy leads, upon a change of variable 8 = 26, to the well-known equa- 
tion of the pendulum, 

(2.7) e -)- r?^ sin e = 0, 

that librates with small oscillations around the laser polarization axis with a 
frequency 

1 f^f^l^'o , l.„c2^'/' 



where the polarizability remains the leading term as soon as the laser frequency 
exceeds 315 cm^^ (for the HCN molecule), whatever the intensity. 

2.3. Control scheme. The question of how to take advantage of the pendular 
mechanism to control molecular alignment can be answered by referring to adiabatic 
vs sudden transport dynamics, starting from an initial isotropic distribution (J = 
M = 0, as prepared, for instance, by laser cooling methods |2]). 

The adiabatic switching of a laser with rise and fall times of 10 ps (same 
order as the molecular rotational period) and an intensity of 10^^ W/cm^ allows 
an adiabatic transport on a single (lowest) pendular state with a good alignment 
during the excitation. Figure |2a) illustrates the alignment dynamics through the 
expectation value of cos^ 0, i.e., 

„2 



(2.9) {cos^0){t)= cos^0\i}{0,t)\sm0d0. 

Jo 

The higher the value, the better the alignment, with <^cos^ 0) — 1/3 corresponding 
to an isotropic distribution. An alignment of (cos^ 0) « 0.64 is achieved during the 
pulse but, as expected, the laser extinction adiabatically transports the pendular 
state back to the isotropic (J — M = 0) state. 
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Figure 2. Alignment dynamics of the HCN molecule submitted 
to (a) a long, adiabatic laser pulse of duration 10 ps and intensity 
10^2 W/cm^; (b) a short, sudden pulse of 1.7 ps duration and 
of intensity 10^^ W/cm^. The pulse envelopes are schematically 
represented by the shaded areas. 

To achieve post-pulse alignment, one may proceed with sudden excitations. 
Figure |2Ib) shows the effect of a sudden pulse of 1.7 ps duration, which excites 
several pendular states with a specific distribution (for a field frequency resonant 
with HCN internal vibrations). The particular J distribution reached is such that 
alignment remains even after the pulse is off, (cos'^ 0) (t) showing small amplitude 
oscillations around a value of 0.5. 

3. Two-color laser control of orientation (Parity breaking scenarios) 

As compared to alignment, orientation, which imposes a given direction, is a 
more challenging goal that requires symmetry breaking in the forward and backward 
directions. Two such scenarios are considered in this work. The first concerns the 
use of two-color pulses, that is a combination of a fundamental frequency to and its 
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Figure 3. Orientation dynamics of the HCN molecule submitted 
to a laser pulse (1.7 ps duration) combining fields of frequency 
Lu and 2w, the latter being in resonance with the v = 1 
vibrational transition of the C-H stretch. 



second harmonic 2oj with a resulting field 

(3.1) £{t) =£o{t)[cosLot + -fcos{2ujt + (p)], 

where maximum asymmetry is reached for 7 ~ 0.5 and = 0. But the asymmetry 
we are referring to is not in the time dependence of £{t), which in any case is not 
enough to produce orientation as, in average over an optical period of the field, 
the positive and negative contributions are canceling each other. We are rather 
referring to parity breaking in the rotational quantum numbers J. As opposite to a 
monochromatic field that excites J's of the same parity, two-color excitation has the 
potentiality to mix J's of different parity. Actually, by taking 2uj in resonance with 
the u = ^ 1 vibrational transition of HCN, the absorption of a single photon 
2uj excites odd J's through the permanent dipole interaction fj,o£, whereas the 
absorption of two photons w excites even J's through the polarizability interaction 
a£^ (starting from an initial state J = M = 0) |Sj. It is precisely the mixture of odd 
and even J's that produces, by the superposition of associated spherical harmonics, 
the asymmetrical angular distributions that are looked for. The result is shown in 
figure|2ta) as the expectation value of cos9, which is a measure of orientation (the 
higher the absolute value, the better the orientation). We see that orientation is 
actually achieved at the extinction of the laser pulse and more efliciently a short 
time after the pulse is off, with (cos0) reaching 0.3. 

Orientation can also be achieved by an adiabatic transport from an initial 
isotropic state to an asymmetric combination of pendular states by a two-color 
excitation scheme |16L I17| . Contrary to the result shown in figure 13 such an ori- 
entation can only be maintained while the laser is on, similarly to what is observed 
for alignment in figure El^a). 
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Figure 4. Electric field of a half-cycle pulse as a function of time 
with its Fourier transform given in the inset. (Adapted from refer- 
ence j28( . taken from C. M. Dion, A. Keller, and O. Atabek, Eur. 
Phys. J. D 14 (2001), 249-255, copyright (2001) by EDP Sciences.) 

4. Half-cycle pulses and the kick mechanism 

The second symmetry breaking scenario addresses the temporal shape of the 
electric field. The most asymmetrical situation that is experimentally reachable is 
the use of half-cycle pulses (HCP's) Such picosecond, far-IR electromagnetic 
pulses are generated when illuminating a GaAs wafer, in the presence of a pulsed 
electric field applied across the surface, by a Ti-Sapphire laser |28| . Such pulses 
have three interesting characteristics (as illustrated on figure 0J: 

(i) their intensities remain less than 10^ W/cm^, therefore present no risk of 
ionization damage on typical diatomic molecules; 

(ii) their relatively broad Fourier transform (^ 30 cm^^) shows the presence 
of non-negligible field components that can induce rotational excitations; 

(iii) they are very sudden with respect to rotation (their positive component 
lasts typically 1 ps) and very asymmetric (the ratio of the positive to the 
negative maximum is more than 10). 

A direct consequence is that such unipolar pulses can impart a "kick" to the mol- 
ecule that effectively orients it in the direction of the positive field component, the 
transfer of angular momentum taking place on the timescale of the short pulse . 
The long-duration negative tail of the HCP is adiabatic enough not to induce fur- 
ther noticeable changes in the dynamics. Actually, the kick seems the most efficient 
orientation mechanism up to date. The results are illustrated on the LiCl molecule 
in terms of the dynamics of the average value of cos 9, in figure [S] 

As in the two-color excitation scenario, the orientation remains modest during 
the positive part of the pulse, but it develops later and reaches the value (cos 9) ~ 
0.5 over a large time period of ~ 4 ps. The kick mechanism by itself can be 
evidenced by working out a sudden impact model Such a model is based on 

the short duration {tp = 1 ps) of the pulse as compared to the molecular rotational 



8 



OSMAN ATABEK AND CLAUDE M. DION 





0.8 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 , 1 1 1 1 , 1 1 1 1 




0.6 






0.4 




A 
CD 


0.2 




<cos 


0.0 
-0.2 






-0.4 






-0.6 






-0.8 






5 10 15 20 25 



time (ps) 

Figure 5. Orientation dynamics of the LiCl molecule submitted 
to the half-cycle pulse of figure 0] (solid line) and using a sudden- 
impact approximation (dotted line). (Taken from C. M. Dion, 
A. Keller, and O. Atabek, Eur. Phys. J. D 14 (2001), 249-255, 
copyright (2001) by EDP Sciences.) 



period (T^ot — 24 ps). The calculation consists in freezing the rotational dynamics 
during the time the positive part of the HCP is acting. At the lowest order, the 
approximation is valid as far as 

(4.1) :^^^o.l3J(J+l)«l, 

where J is the rotational kinetic operator and B the rotational constant. The ex- 
cellent agreement between this sudden impact approximation and the exact time- 
dependent calculation, with a root-mean-squared difference less than 0.03, elo- 
quently advocates for the kick mechanism. 



5. Optimal control scenarios 

5.1. General frame. A completely different approach in its philosophy, with- 
out a priori reference to any dynamical mechanism, is optimal control. Such a 
scheme presents a general structure resting on three steps: 

(i) an evaluation function, {cos 9) (t), involving the numerical solution of the 
time-dependent Schrodinger equation (using split operator |3l 1131 114| or 
basis expansion techniques) . For a laser-driven rigid rotor this amounts 
to solving 

(5.1) ih^^^{e,v,t)^H{t)^{e,>p,t) 

with 

(5.2) H{t) ^ BJ"^ - /iof (i) cos 9 - [{a^ - a^^) cos^ + aA_] 



MOLECULAR ALIGNMENT AND ORIENTATION 



9 



in order to calculate 

/'27r ('IT 

(5.3) {cose){t)= / \ij {9, (p; t) f cos sine d9d(p 

Jo Jo 

(ii) a set of parameters defining the temporal shape of the laser pulse, given 
here as a superposition of individual sine-square functions with interme- 
diate plateau values, 

N 

(5.4) £(t) = ^£„(t)sin(w„t-|-0„) 



n=l 



with 



(5.5) 



£n{t) = < 
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(iii) a target j, i.e., an optimization criterion. 

The optimization aims at the obtainment of the laser parameters (intensity, 
frequency, temporal shape defined by rise and fall times, together with possible 
plateau durations) that optimally satisfy the target j. The optimization proce- 
dure we use is based on genetic algorithms [D 122) operating in a 7N dimensional 
parameter space (i.e., 7 parameters for each of the N individual pulses). It is 
worthwhile noting that a given criterion puts the emphasis on a specific feature 
of the orientation, such that different criteria may lead to quite different results. 
This is basically related to the fact that orientation can by no way be achieved by 
monitoring a single quantum state or a unique and pre-determined superposition 
of states. We will analyze in some detail the consequences of the choice of specific 
criteria. But before doing this, let us proceed to a first optimization retaining an 
intuitive, simple criterion that consists in the maximization of {cos 9) for a single 
time tf when the laser field is off, 

(5.6) j = \{cos9){ti)\. 

5.2. The kicked molecule. The optimal electric field which is obtained and 
illustrated in figure El^a) turns out to be one of the most enlightening results of 
this study. A combination of three individual pulses (i.e., 21 parameters to be 
optimized) leads to a very sudden and asymmetric pulse, followed by a null plateau 
and some rather symmetrical oscillations. The molecule (HCN taken as a rigid 
rotor, in this example) orients through the already discussed kick mechanism due 
to the fast angular momentum transfer by the unipolar pulse at about t = 0.3 ps 
and during the null plateau lasting 1 ps. The symmetrical oscillations of this field 
between 1 and 1.7 ps do not seem to play a major part in this dynamics. Here 
again, the kick mechanism can be evidenced by a sudden impact approximation 
which closely follows the {cos 9) dynamics at least until 1.5 ps [see figure Elb)]. 
It is important to note that an accuracy of 4 digits in the genetic algorithm is 
achieved for the intensities, frequencies, and absolute phase differences in order to 
produce the particular shape of the electric field of figure E^a) , responsible of the 
subsequent kick mechanism. Figured accounts for the long-time dynamics (i.e., for 
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Figure 6. (a) Laser field resulting from the optimization us- 
ing criterion j = \{cos9) {tf)\, for the HCN molecule in a rigid- 
rotor model, (b) Orientation dynamics obtained using this field 
(solid line) and using a sudden-impact approximation (dashed 
line). (Taken from C. M. Dion, A. Ben Haj-Yedder, E. Cances, 
C. Lc Bris, A. Keller, and O. Atabck, Phys. Rev. A 65 (2002), 
063408, copyright (2002) by the American Physical Society.) 



full rotational period). Maximum orientation is reached at about t « 10 ps and 
lasts for about 3 ps. It is interesting to note that the result of the sudden impact 
model based on the first kick (at ^0.5 ps) can further be improved by introducing 
the effect of the oscillatory field acting after the null plateau as a second kick (that 
is again freezing the rotational motion during the time over which the molecule 
experiences this second interaction). 

An analysis of the optimized laser pulse obtained can be carried out using the 
short-time Fourier transform of the electric field, given by 



(5.7) 



/+00 
-oo 



£{T)U{T-t)e-'^-dT, 
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Figure 7. Same orientation dynamics as figm-elHtb), but given for 
an entire rotational period of the HCN molecule (Trot = 11.45 ps) 
and with an additional curve (dotted line) corresponding to the 
sudden-impact approximation with two kicks (see text). (Taken 
from C. M. Dion, A. Ben Haj-Yedder, E. Cances, C. Le Bris, 
A. Keller, and O. Atabek, Phys. Rev. A 65 (2002), 063408, copy- 
right (2002) by the American Physical Society.) 
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Figure 8. Short-time Fourier transform [equation 1)5. 7|l ] of the 
field given in figure EI a). 



with /w a windowing function (we have taken a Tukey-Hanning window |21| of 
total width 0.29 ps). The result is shown in figure |S1 as a 2D representation of the 
electric field amplitude as a function of w and t. The asymmetric feature present in 
the field at around ~ 0.15 ps turns out to be the superposition of two components 
at frequencies uj and 2a;. This corresponds to the most asymmetrical electric field 
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which, when truncated at a time that is not an integer of the optical period, leads 
to a non-zero field integral imparting the kick to the molecule. This is followed later 
on by a component centered around a single frequency which also leads to a second, 
smaller kick. Such a result has a broader applicability: the same pulse applied to 
a different molecule (e.g., LiF) produces similar orientation dynamics [9]. 



5.3. The optimization criteria. We now return to a thorough analysis of 
the consequences of choosing different criteria in the optimization scheme. Ideally, 
the best orientation is the one which is the most efhcient (i.e., with a value of 
(cos 6) close to one) and which lasts forever. The relative merits of the different 
criteria can thus be quantitatively measured in terms of the efficiency and the 
duration of the orientation. All calculations that are presented hereafter apply to 
the LiF molecule and are conducted within a parameter sampling space limited to 
frequencies comprised in between 500 and 4000 cm^^, to individual pulses extinction 
times less than ti taken as one tenth of a rotational period (as we are seeking post- 
pulse orientation resulting from a sudden excitation, which also helps in keeping 
calculations within acceptable CPU times) and to laser intensities not exceeding 
3 X 10^^ W/cm^ (to avoid risk of ionization damage, having in mind that 1.1 x 
10^* W/cm^ is the predicted first ionization threshold for LiF). 

The criteria can roughly be classified into two groups: simple, when they em- 
phasize either efficiency or duration; hybrid, when they realize a compromise be- 
tween efficiency and duration. Aiming at the best orientation efficiency, a more 
flexible criterion then the one that has been used so far [equation 1)5. is 

(5.8) ji= max | (cos 6*) (i)| , 

i6[t/,t/+Trot] 

the maximum being now taken at any time within the interval [i/, + ri.ot] (taking 
into account the periodical behavior of the orientation dynamics with respect to 
the rotational period Trot)- The result is displayed in figure |S1 The electric field 
thus obtained is built up from two individual pulses of comparable intensities, with 
frequencies in a ratio of 2 responsible for the double wiggles in the amplitude, as 
illustrated in figure [TT)I Excellent orientation {{cosff) w —0.8) is achieved with a 
time delay of about 6 ps after the pulse is off. This molecular response time to 
the laser interaction is the time needed for the specific phase interferences to occur 
in the rotational wave packet. Such occurrences of the field- free dynamics are 
periodical (molecular rotational periodicity Trot), give rise to the revival structures 
discussed in the literature j25| . and are common to all results presented hereafter. 

The criterion ji fails in producing long-lasting orientation (the time interval 
over which | (cos 0) \ remains larger than 70% of its maximum value does not exceed 
0.37 ps). Another criterion is built to precisely maximize the orientation duration 

(5.9) 32 = max 

te[tf,tf+T,^t] Trot 

where t is defined as the time over which an orientation exceeding is kept, 

(5.10) iL<|(cos0)(t)|<ii. 
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Figure 9. (a) Laser field resulting from the optimization using 
criterion ji — majiteltf.tf+Trat] I(cos0)(t)|, for the LiF molecule 
in a rigid-rotor model, (b) Corresponding orientation dynam- 
ics. (Adapted from A. Ben Haj-Yedder, A. Auger, C. M. Dion, 
E. Cances, A. Keller, C. Le Bris, and O. Atabek, Phys. Rev. A 
66 (2002), 063401, copyright (2002) by the Americal Physical So- 
ciety.) 



The optimized electric field shows two intense, sudden, high frequency and well 
separated components (see figure , inducing a smooth behavior for (cos 6) (t) as 
expected for large values of t, but unfortunately the orientation is very inefficient. 

The weakness of these simple criteria may be overcome by building physically 
more sound hybrid criteria that combine the advantages of the previous require- 
ments for a better compromise between efficiency and duration. Our best result is 
obtained with a criterion that maximizes a functional of {cos 9)^ (t) over an entire 
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Figure 10. Short-time Fourier transform [equation (|5.7|) ] of the 
field given in figure a). 



rotational period, 



(5.12) C{t) = 



(5.11) j3 = niax-^ / C'it)dt, 

J rot Jtf 

where C(t) is tailored as to put the emphasis on time intervals where {cos 9) is larger 
than some fixed value (0.4 in our case) through an appropriate weighting factor (0.1 
in our case): 

0.1 (cos 6*) {t) if (cos 6*) {t) < 0.4 
{cos 6) (t) elsewhere 

Figure 1121 displays the optimized field built up from two sudden, intense pulses, 
with very close frequencies (742 cm^^ and 808 cm^^) giving a beat-like structure. 
It is worthwhile noting that the resulting field appears as two beats of similar 
frequency ~ 780 cm~^ (see figure [T^ . but we failed in analyzing them in terms 
of a double kick mechanism (the sudden impact model not providing the expected 
orientation dynamics). Finally, in terms of quantitative measurements, this field 
achieves orientation which lasts for 0.7 ps with a |(cos6')| larger than 0.5, and is 
one of the best results of the literature. 

5.4. Robustness with respect to temperature. All calculations so far 
presented concern an initially rotationless molecule, J = M = 0, M being the 
quantum number labelling the projection of the rotational angular momentum on 
the laser polarization vector taken as the laboratory reference axis. With respect 
to experimental conditions, such a situation concerns a molecule at a temperature 
of T = K. Starting from such a state (prepared by laser cooling methods, for 
instance |2]), linearly polarized photon absorption only populate higher J levels; 
M, being a good quantum number (labelling an azimuthal angular motion, well 
separated due to cylindrical symmetry), remains unchanged (i.e., zero). Increasing 
temperature, on the other hand, populates higher J levels with all their M compo- 
nents, \M\ — Q,l, . . . , J, within a Boltzmann distribution. For a given component. 
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Figure 11. (a) Laser field resulting from the optimization us- 
ing criterion ji — max^g^j^. jj^+t^^j r/Trot, for the LiF molecule 
in a rigid-rotor model, (b) Corresponding orientation dynam- 
ics. (Adapted from A. Ben Haj-Yedder, A. Auger, C. M. Dion, 
E. Cances, A. Keller, C. Le Bris, and O. Atabek, Phys. Rev. A 
66 (2002), 063401, copyright (2002) by the Americal Physical So- 
ciety.) 



a vectorial representation of the field-driven classical rigid rotor leads to 

\M\ 



(5.13) 



smf = 



the total angular momentum vector J being orthogonal to the rigid rotor axis R and 
d being the polar angle positioning R with respect to the laser polarization vector. 
Large values of M would thus prevent alignment and orientation, Q reaching a value 
as high as 7r/2 for \M\ = J, as illustrated in figurelTH Using linearly polarized laser 
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Figure 12. (a) Laser field resulting from the optimization us- 
ing criterion js = max (l/Ti-ot) J*^^ dt, for the LiF mole- 
cule in a rigid-rotor model, (b) Corresponding orientation dynam- 
ics. (Adapted from A. Ben Haj-Yedder, A. Auger, C. M. Dion, 
E. Cances, A. Keller, C. Le Bris, and O. Atabek, Phys. Rev. A 
66 (2002), 063401, copyright (2002) by the Americal Physical So- 
ciety.) 



sources, no reduction of M could be achieved, showing that orientation is no more 
a mathematically controllable process at non-zero temperature. Because they are 
breaking the overall cylindrical symmetry, elliptically polarized lasers are possible 
tools for modifying the values of M, in an optimally controlled way to lead again to 
alignment and orientation. Work in this direction is now undertaken in our group. 

It remains however that the robustness of our optimal control scheme can be 
checked at least for some low temperatures. Instead of our previous evaluation 
function H5.3|l . we have to refer to a temperature dependent Boltzmann averaging 
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Figure 13. Short-time Fourier transform [equation (|5.7|) ] of the 
field given in figure ^fa). 
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Figure 14. Schematic representation of the total angular momen- 
tum vector J and its projection M on the laser polarization axis 
£, for a (diatomic) linear molecule. By definition, J is perpendic- 
ular to the internuclear axis, therefore if |M| = J the molecule is 
perpendicular to the laser polarization axis. 

of cos 9, 

(5.14) ((cos 61)) (t) ^ Q-^ exp 

,/ 

where (cos d) j (t) results from the dynamics of an individual initial state J, M 
included in a summation over J and M with its appropriate exponential weighting 
factor, and a partition function 

Jmax 

(5.15) Q^Yl (2-/+ l)exp 

J 
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for normalization. 

Figure [T5l presents the results for the LiF molecule at T = 5 K. Here again, 
simple (ji) and hybrid (js-like) criteria are referred to. The dotted line (panel a) 
represents the orientation dynamics at 5 K [equation H5.14|l ] but with the electric 
field of figure IHt a) optimized for T = K. This is, precisely, to show how orienta- 
tion effects are rapidly washed out when incorporating higher initial M's. A field 
optimized for T = 5 K restitutes again the oscillatory structures of {{cos 9)) (t), 
with however reduced amplitudes (solid line of panel a). A more satisfactory re- 
sult requires the introduction of an additional individual pulse bringing 21 free 
parameters in the optimization scheme. The result in given in figure |^b) with an 
orientation higher than {{cos9)) = 0.27 kept over a time of about 0.25 ps. The re- 
sult of the hybrid criterion j3 with a 21-free-parameters optimization is illustrated 
in figure ini[c). A much longer time duration (~ 0.4 ps) is achieved for about the 
same efficiency. 



6. Conclusion 

The complementarity between basic mechanisms and optimal control schemes 
can be summarized either by referring to the way of implementing mechanisms in 
control, or reversely of identifying mechanisms revealed from control. 

The first strategy rests on three observations: 

(i) there is no single solution arising from optimal control scenarios (not only 
different criteria, but also different sampling spaces for parameters would 
lead to different results); 

(ii) a careful study of the laser-induced dynamics can help to the identification 
(or guess) of some basic mechanisms leading to the desired observable; 

(iii) by appropriately building the targets and delineating the parameter sam- 
pling space, we can help the optimal control scheme to take advantage of 
the mechanisms. 

The second strategy consists in learning a dynamical mechanism presenting 
some robustness from the optimal control results. This mechanism may be used 
directly for reaching the desired observable (strategy 2) or further improved by 
implementing it again in a control scenario with some additional flexibility (strategy 
1). _ 

The kick mechanism and the optimal control carried out using sudden pulses 
is an example that illustrates these two routes. The mechanism has actually been 
depicted using half-cycle pulses and a sudden impact model. The parameter sam- 
pling space for the optimal control, although not constrained, has been limited to 
sudden pulses that are precisely compatible with the kick mechanism. This control 
scenario has finally revealed a double-kick mechanism which can serve as a starting 
point of a train of kicks, that we can merely guess to provide even better orienta- 
tion. Finally, this mechanism, together with some additional free parameters such 
as the time delay between the successive kicks, can be re-injected (re-introduced) 
in an optimal control scheme (strategy 1) aiming at improving its efficiency. We 
are presently pursuing some investigations in this direction. 
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Figure 15. Effect of temperature on the orientation dynamics of 
tlie LiF molecule, (a) Dotted line: using the laser field obtained 
from a T = K optimization with criterion ji [see figure|^a)], but 
applied to an initial distribution of molecules at T 5 K. Solid 
line: using the laser field obtained from a T = 5 K optimization 
with criterion ji. (b) Same as (a), but using a superposition of 3 
laser pulses, (c) Same as (b), but using criterion js. (Adapted from 
A. Ben Haj-Yedder, A. Auger, C. M. Dion, E. Cances, A. Keller, 
C. Le Bris, and O. Atabek, Phys. Rev. A 66 (2002), 063401, 
copyright (2002) by the Americal Physical Society.) 
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